Josephson effect as a measure of quantum fluctuations in the cuprates 
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Quantitative analysis of the Josephson effect is shown to provide direct information about phase 
fluctuations in the superconducting banks. Applying the analysis to the cuprates, substantial quan- 
tum fluctuations between d-wave and s-wave pairing are found at low temperatures. A phenomeno- 
logical model of such fluctuations is introduced and solved in a self-consistent single-site approx- 
imation. The cuprate parameters are argued to lie close to an exotic superconductor phase with 
condensation of electron quadruplets. 



More than a decade ago, it has been suggested that 
the Josephson effect can be used as a direct probe of the 
angular dependence of the order parameter phase in the 
cuprates Since then, several types of phase sensitive 
experiments have been performed, and an overwhelming 
majority of them confirmes the d-wave symmetry of pair- 
ing in the cuprates (for a review, see 0). 

In absence of time reversal symmetry breaking, the 
superconducting current / as a function of the phase dif- 
ference (j) across the Josephson junction reads — 
J2n sin ncj) , and in the usually apphcable tunnel 
limit only the first harmonics Ii needs to be considered. 
In typical phase sensitive experiments , one tests only 
the dependence of the sign of /i on the junction geome- 
try. In a recent series of experiments the function I{4>), 
which contains a more complete information about the 
equilibrium superconducting properties of a Josephson 
junction, has been measured for junctions involving the 
cuprates (for a review, see 0). 

The present work is motivated by the special features 
of /((/>) found in the cuprates, in particular by the small 
first harmonics and the large second harmonics of I{4>). 
The outline of this Letter is as follows. First we show 
that a quantitative analysis of the current-phase relation 
I{4>) provides information about the phase fluctuations 
in the superconducting banks 0, 0| . Further we intro- 
duce a phenomenological model of phase fluctuations in 
the cuprates and solve it within the simplest nontrivial 
self-consistent approximation. Finally we discuss the mi- 
croscopic justification of the phenomenological model. 

Let us start by considering a Josephson junction de- 
scribed, in absence of phase fluctuations in the banks, by 
the harmonics /„. Now, let us switch on phase fluctua- 
tions and denote the fluctuating part of the phase in the 
bank i = 1,2 as ipi. In the tunnel limit the phase fluctu- 
ations in different banks can be considered independent 
of each other. Thus at time scales longer than the typi- 
cal fluctuation time and/or for junctions larger than the 
fluctuation correlation length, /(0) renormalizes to 

H4>) = ^/,i(sinn(0 + ipi~ 1P2)) = sinn0, 

n n 

where the renormalized n-th harmonic reads /«//„ = 
{cos mpi) (cos nip 2). In grain boundary junctions with a 




FIG. 1: (a) Average phase field configuration (ipi) (circles 
at the dual lattice sites) in a d-wave superconductor. Note 
the two sublattice structure of the dual lattice, (b) A typical 
snapshot of the fluctuating phase field (pi — {(pi) + Si. A 
lattice site 7? (or the corresponding lattice plaquette) is said 
to support local d, —d, s, and —s pairing, if phases at the 
bonds {R,R + x) and {R,R + y) are (+,-), (-,+), (+,+), 
and (—,—), respectively. 

45° misorientation between the cuprates, the first two 
harmonics are experimentally accessible. A detailed anal- 
ysis of the data 4, 5] suggests that 

a = (cost/?) « 0.3, /? = (cos2(p) « 1. (1) 

We would like to point out that the value of a was deter- 
mined from the Josephson product IqRn of a large set 
of cuprate/cuprate and of cuprate/low- Tc junctions (for 
reviews, see 0, 01) which are renormalized by and 
a, respectively. The values of a determined using both 
types of junctions are in semiquantitative agreement with 
each other and as such appear to be quite reliable 0, . 
On the other hand, the value of f3 is subject to a larger 
uncertainty, being based on a much smaller set of exper- 
imental data 0. However, the (anomalous) inequality 
f3 > a seems to be quite reliable. 

Now we proceed by introducing a minimal model of 
phase fluctuations consistent with Eq. (Q. Our model 
is motivated by the strong-coupling RVB picture of the 
cuprates which views the cuprates as a liquid of predom- 
inantly local singlets While the pairs are formed at 
a high energy scale ~ J, they acquire phase coherence 
only at a much lower scale. Therefore we assign a fluc- 
tuating phase field ipi to each bond i of the square Cu 
lattice. Note that (pi is the phase of a Cooper pair, i.e. 
of a charge 2e object. 

Let us first consider a toy model in which the fluc- 
tuating part 5i of the phase at each link is either 
or ±7r, with probabilities Pi = (1 -f Q:)/2 « 0.65 and 
P2 = (1 ~ Q^)/2 ~ 0.35, respectively. One finds readily 
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that Eq. is satisfied in such a model. On the other 
hand, Fig. ^ shows a snapshot of a typical ipi configu- 
ration for such fluctuations. Note that the probabilities 
that a lattice point supports local d, —d, s, and —s pair- 
ing symmetry are Pi « 0.42, P| « 0.12, P1P2 « 0.23, 
and P1P2 ~ 0.23, respectively. In other words, the nu- 
merical values in Eq. |^ suggest the presence of strong 
quantum fluctuations from d-wave to s-wave pairing 

The simplest model of a superconducting Cu02 plane 
which allows for phase fluctuations, takes into account 
the compact nature of the phase field, and assumes dom- 
inant d-wave and subdominant s-wave pairing (with pair- 
ing strengths V + W and V — W, respectively), reads 

H = -m-'Y.9'/9v' + Y.^■^iV'^,VJ)■ (2) 

i (ij) 

The first term in Eq. ^ which allows for quantum fluc- 
tuations of phase can be thought of as a local charging 
energy at site i. Note that from now on, instead of bonds 
we talk about the sites of a (dual) lattice. The second 
term where (i, j) denotes a pair of nearest neighbor sites 
can be viewed as a Josephson coupling between the sites i 
and j, Ej{(pi,ipj) = —Vcos{2ipi — 2ipj) + Wcos{Lpi — (pj). 
Note that the sign of the W term is different from the 
standard convention, because of the d-wave symmetry of 
pairing in the cupratcs. 

Let us point out that our model Eq. ||2Jl includes the 
quantum XY model Q as a special limit V — 0; we will 
show that the most spectacular effects occur at V ^ W. 
Note also that we haven't considered dissipative effects in 
Eq. (0). Since the cuprates are d-wave superconductors 
with low- lying fermionic excitations , this is a nontriv- 
ial assumption which needs to be studied in future work. 

In what follows we study the model Eq. ^ by con- 
structing a variational solution of the form 

■0(<y5l,<^2, ■ ■ •) = niVi(</'»), (3) 

where the product runs over all lattice sites and the local 
functions ipi are 2tt periodic. The Schrodinger equation 
for the lattice point i reads 

[-m-i^W + U,iip)]Mip) = e,Mv), (4) 

with a self-consistent potential 

Ui{ip) = ^ [-~V{cOa2(p)i^r C0s2ip + W {cos (p)i+r COS ip] , 
r 

where the sum runs over the four nearest-neighbor sites 
i + T of the studied site and {f{<p))i = Jdipf{ip)\ilJi{ip)\'^. 
Note that the single-site approximation Eq. (O is similar 
in spirit to the self-consistent harmonic approximation 
, but goes beyond it by treating the local problem ex- 
actly, which is instrumental in the present context. 

Now we specialize to the case of a two-sublattice solu- 
tion with '4>i(ip) = iPa{^) and = tpsif) for i in sub- 
lattice A and B, respectively. Note that this corresponds 
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FIG. 2: Solid line: Potential energy entering Eq. |^ (in arbi- 
trary units). Result of a self-consistent calculation for « = 1.1 
and w — 0.09, leading to a « 0.22 and b « 2.0. Dashed line: 
Square of the ground-state wavefunction |i/)((^)|^ for the same 
parameters. 

to a translationally invariant solution on the original lat- 
tice. Assuming furthermore a d-wave solution, iPb{^) = 
^JAif + tt), we obtain {cosn(p)B = {— 1)"" (cos nip) a- 
The self-consistent Hamiltonian for sublattice A there- 
fore reads 

H = —d'^ /dip'^ — a cos If — bcos2tp, (5) 

where we have measured the energy in units of m^^. In- 
troducing dimensionless interaction parameters v — mV 
and w — mW, the dimensionless self-consistent poten- 
tials read 

a = Awa, b = 4:v(3, (6) 

where a and /? are defined in Eq. in which the mean 
values (...) are to be calculated with respect to the solu- 
tion of Eq. (O . The total energy per lattice site £ differs 
from the lowest eigenvalue e of the Hamiltonian Eq. jSjl: 

e = e + {aa + bp)/2. (7) 

In Fig.|2]we plot the potential energy entering Eq. Q 
in the limit b ^ a relevant to the cuprates. The classical 
solutions at if = and ip = ±7r are nearly degenerate, 
their splitting being 2a. The ground-state wavefunction 
Tpiip) is localized predominantly in the vicinity of ip = 0, 
but tunneling across the finite barrier leads to a second 
peak of ipif) in the vicinity of tp — ±7r. In complete 
analogy with the toy model Fig. ^ this tunneling leads 
to quantum fluctuations between c?-wave and s-wave pair- 
ing, which are thus seen to be essentially local. 

Self-consistent solution to Eas. H5l()ll|l shows that three 
phases exist in the w > part of the phase diagram of 
the model Eq. (0), see Fig.|21 (i) a conventional d-wave 
superconducting phase with a > and /3 > realized for 
sufficiently large v,w; (ii) a quantum disordered phase 
with a ~ /3 — at small v,w; and (iii) an exotic su- 
perconductor phase with a — and /3 > appearing at 
large v and small w. 
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FIG. 3: Phase diagram of the model Eq. ((2j. We hypothe- 
size that superconducting cuprates track the hne between the 
points OD and UD, corresponding to overdoped and under- 
doped materials. The dotted lines show possible trajectories 
under further reducing the hole density. 

Extending the phase diagram for w > to negative val- 
ues of w, one verifies easily that the role of cZ-wave and 
s-wave pairing is interchanged and therefore the phase 
diagram for w < is a mirror image of that for w > 0. 
This implies that within the quantum model Eq. |(2J), 
upon changing the sign of W the superconducting state 
switches from d-wave to s-wave symmetry by crossing 
one of the two nontrivial phases in Fig. |31 Which of 
the two possible intermediate states is realized depends 
on the strength of quantum fluctuations. In the limit of 
strong fluctuations (small v), an intermediate quantum 
disordered state is realized [phase (ii)], which is presum- 
ably insulating. If the strength of quantum fluctuations 
decreases, the exotic phase (iii) separates the two conven- 
tional pairing states. In the classical limit m — > cxo the 
width of phase (iii) vanishes. In what follows we argue 
that although {cos(p) = in phase (iii), it does represent 
a superconductor, albeit an exotic one, since it can be 
thought of as a condensate of pairs of Cooper pairs, or 
quadruplets of electrons . Note that phase (iii) bears 
some similarity to the fractionalized superconductor dis- 
cussed in ^jj. 

Let us show now that phase (iii) is described 
by a macroscopic Ginzburg-Landau type wavefunction 
^GLix,y) with a nonvanishing phase rigidity. To this 
end consider the many body wavefunction Eq. Q with 
the following choice of local wavefunctions, 

'^P^+vif) = tpi{V + '^ + X2), (8) 

where x and y are unit vectors of the lattice in the x 
and y directions, respectively. Eq. © corresponds to a 
macroscopic wavefunction with an internal d-wave sym- 
metry and a spatially modulated phase '^GL{x,y) oc 
exp —i(xix + X'iy) / do, where do is the lattice constant 
of the dual lattice. The Schrodinger equation for site i 
can still be written in the form of Eq. 10} with a selfcon- 
sistent potential satisfying [/^^^((p) — Ui{(f -\- t: + xi) 



FIG. 4: Order parameters a and /3 along two cuts of the 
phase diagram Fig.|3 v — 1.5 (solid lines) and v = 2 (dashed 
lines) . The inset shows the total energy per lattice site £ as a 
function of v along the line «; = 0. 

and Ui+y{ip) = Ui{(p + TT + X2)- The periodicity of 
the local wavefunctions together with Eq. © guaran- 
tees that all local energies Si are equal. For a gen- 
eral lattice point i the potential Ui{ip) contains terms 
proportional to both cos nip and sin nip with n = 1,2. 
However, we assume that there exists one special lattice 
point for which the terms propoportional to sin nip van- 
ish. For this lattice point the Hamiltonian can be writ- 
ten in the form of Eq. jSJ. The set of self-consistent 
equations is closed by Eq. in which the effective 
interaction strengths Wcft = u>(cosxi + cosx2)/2 and 
VcS = u(cos2xi + cos2x2)/2 have to be used instead 
of w and v, respectively. Since the total energy per lat- 
tice site is still given by Eq. |7J, the energy difference 
between a modulated and a uniform macroscopic state 
can be written in a Ginzburg-Landau type form 

ds dS 

&£ = e(woff , Woff) - £{w, v) w -^5w - —5v, (9) 

ow ov 

where 5w — w — Wcs ~ w{x\ + xi)/^ and 5v = v — v^s ~ 
d(Xi + xi) and the approximate equations are valid for 
slowly varying macroscopic wavefunctions, x1jX2 -C 1. 
In the conventional superconductor phase [phase (i)] we 
expect in a generic point of the phase diagram de/dw < 
and de/dv < 0. From Eq. © it then follows that in the 
long- wavelength limit Se oc (xi + xDi as it should be in 
a superconductor with a finite phase stiffness. On the 
other hand, in phase (iii) we have de/dw = 0, but the 
inset to Fig. ^ shows explicitly that de/dv < 0, which is 
enough to guarantee a finite phase stiffness in this case 
as well. Thus we have shown that the phase (iii) is an 
exotic superconductor. 

Returning to the physics of cuprates, Eq.l^ implies 
that if we insist on their description in terms of the model 
Eq. (0), then its parameters should be chosen inside the 
phase (i), but close to phase (in). Fig. ^ shows explic- 
itly that in this case it is possible to have a « 0.3 < (3. 
Within our simplistic model it is difficult to have « 1 
at the same time, because this would require small quan- 
tum fiuctuations (large v). However, for f 3> 1 the differ- 
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FIG. 5: (a) Superconducting transition temperature of the 
t-t' Hubbard model in the d-wave {Ted, squares) and s-wave 
{Tcs, circles) symmetry sectors for t' /t — 0.45 and U/t — 6 
as a function of electron density. Also shown is the transi- 
tion temperature Tcxy in the leading subdominant symmetry 
sector {dxy-wave, triangles). Calculated following |lal on a 
lattice 1024x1024. (b) Coupling constants V and W (in ar- 
bitrary units) estimated from Ted and Tea- 

ence W between the d-wave and s-wave pairing strengths 
needs to be implausibly small in order to stabilize the ex- 
otic superconductor phase. 

Next we address the question whether it is reasonable 
to assume that, within a realistic model of the cuprates, 
d-wave and s-wave pairing are competitive in energy. 
Lead by the hypothesis that the symmetry of pairing is 
correctly determined already at weak coupling we 
have studied the superconducting phase diagram of the 
t-t' Hubbard model for t' /t = 0.45 and a moderate in- 
teraction strength U/t = 6 as a function of the electron 
density p in the vicinity of p — 0.7. This choice of param- 
eters is close to the canonical values for the cuprates and 
it is inspired by the results of [l^ , where a d-wave/s-wave 
pairing transition was found in this region. 

In Fig.|Slwe plot the superconducting transition tem- 
peratures determined by the variational method for a 
Kohn-Luttinger superconductor [l^ in the two dominant 
pairing symmetries (d-wave and s-wave). Note that for 
extremely overdoped cuprates with p < Pc~ 0.7, s-wave 
pairing is the leading instability. It is also shown that 
close to pc Tc in the leading subdominant pairing sym- 
metry sector ((ia;j,-wave) is much smaller than Ted, Tcs- 

According to Eq. (|2Jl, the condensation energy per site 
in absence of phase fluctuations is 2{V+W) and 2{V—W) 
in the d-wave and s-wave pairing sectors, respectively. At 
weak coupling, these condensation energies are propor- 
tional to T^^ and T^^. By comparing these expressions, 
we obtain a rough estimate of the coupling constants V 
and W. The results are shown in Fig. \E\ Note that 
with hole doping x = 1 — p decreasing from Xc ~ 0.3 
towards experimentally accessible values, both W and V 
and their ratio W/V increase. On the other hand, also 
the fluctuations are expected to grow with underdoping, 
i.e. m^Oasx^O^ This brings us to the 



question about the doping dependence of ?; = mV and 
w = mW. We hypothesize that with reduced doping, the 
effect of TO — > dominates and the cuprates track the line 
in Fig. 131 between the points OD and UD. Under further 
reducing the hole density, the insulating quantum disor- 
dered phase should be reached. This can happen either 
directly or via crossing the exotic superconductor phase, 
as indicated by the dashed lines in Fig. El 

It should be remarked that the question about how 
/(0) depends on the doping level is very difficult to ad- 
dress experimentally. An interesting possibility might be 
to study artificially doped grain boundary junctions [l5| . 

In conclusion, we have found that the two seemingly 
independent experimental observations in the cuprates, 
namely the large first and the small second harmonics 
of the Josephson current, can be explained by a single 
assumption of substantial quantum fiuctuations between 
d-wave and s-wave pairing at low temperatures. We have 
introduced a phenomenological model of such fiuctua- 
tions and solved it in a self-consistent single-site approxi- 
mation. We have argued that the cuprate parameters lie 
close to an exotic superconductor phase with condensa- 
tion of electron quadruplets. Finally, d-wave and s-wave 
pairing were shown to be competitive within the Hub- 
bard model for parameters close to the canonical values 
for the cuprates. 
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